fid1 = fopen('..\..\output\initial_pops.csv');
fid2 = fopen('..\..\output\afterHarvesting.csv');
fid3 = fopen('..\..\output\afterHunting.csv');

A  = fscanf(fid1, '%d\t%d\t%d', [3, inf])';
B  = fscanf(fid2, '%d\t%d\t%d', [3, inf])';
C  = fscanf(fid3, '%d\t%d\t%d', [3, inf])';

rabbitsPopulation = serialize(A(:, 1), B(:, 1), C(:, 1));
wolvesPopulation = serialize(A(:, 2), B(:, 2), C(:, 2));
roesPopulation = serialize(A(:, 3), B(:, 3), C(:, 3));

initX = (0:size(A(:,1))-1) * 3 + 1;
harvestX = (0:size(A(:,1))-1) * 3 + 2;
huntX = (0:size(A(:,1))-1) * 3 + 3;

fclose(fid1);
fclose(fid2);
fclose(fid3);
close all;
figure(1)

%subplot(2,1,1)
hold on;
plot(rabbitsPopulation, 'b'); 
plot(wolvesPopulation, 'g');
plot(roesPopulation, 'r');

scatter(initX, A(:,1), 30, 'o', 'b')  %populacja krolikow w init
scatter(harvestX, B(:,1), 30, 'x', 'b') %populacja krolikow po harvest
scatter(huntX, C(:,1), 30, '+', 'b') %populacja krolikow po hunt 

scatter(initX, A(:,2), 30, 'o', 'g') %populacja wilkow w init
scatter(harvestX, B(:,2), 30, 'x', 'g') %populacja wilkow po harvest
scatter(huntX, C(:,2), 30, '+', 'g') %populacja wilkow po hunt 

scatter(initX, A(:,3), 30, 'o', 'r') %populacja saren w init
scatter(harvestX, B(:,3), 30, 'x', 'r') %populacja saren po harvest
scatter(huntX, C(:,3), 30, '+', 'r') %populacja saren po hunt 

legend('Rabbits', 'Wolves', 'Roes', 'init', 'harvest', 'hunt' );
title('Populacje w czasie (jeden cykl to 3 pkt na osi x)');

run printScores
run printGenetics
